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Abstract 

In this paper we develop general LP and ILP techniques to find an approximate solution 
with improved objective value close to an existing solution. The task of improving an ap- 
proximate solution is closely related to a classical theorem of Cook et al. [7] in the sensitivity 
analysis for LPs and ILPs. This result is often applied in designing robust algorithms for 
online problems. We apply our new techniques to the online bin packing problem, where 
it is allowed to reassign a certain number of items, measured by the migration factor. The 
migration factor is defined by the total size of reassigned items divided by the size of the 
arriving item. We obtain a robust asymptotic fully polynomial time approximation scheme 
(AFPTAS) for the online bin packing problem with migration factor bounded by a polyno- 
mial in ^. This answers an open question stated by Epstein and Levin |10| in the affirmative. 
As a byproduct we prove an approximate variant of the sensitivity theorem by Cook at el. 
[7] for linear programs. 



1 Introduction 

The idea behind robust algorithms is to find solutions of an optimization problem that are not 
only good for a single instance, but also if the instance changes in certain ways. Instances change 
for example due to uncertainty or when new data arrive. With changing parameters and data, 
we have the efi'ort to keep as much parts of the existing solution as possible, since modifying 
a solution is often connected with costs or may even be impossible in practice. Achieving 
robustness especially for linear programming (LP) and integer linear programming (ILP) is 
thus a big concern and a very interesting research area. Looking at worst case scenarios, how 
much do we have to modify a solution if the LP/ILP is changing? There is a result of Cook 
et al. [7] giving an upper bound for ILPs when changing the right hand side of the ILP. Many 
algorithms in the theory of robustness are based on this theorem. 

As a concrete application we consider the classical online bin packing problem, where items 
arrive over time and our objective is to assign these items into as few bins as possible. The 
notion of robustness allows to repack a certain number of already packed items when a new 
item arrives. On the one hand we want to guarantee that we use only a certain number of 
additional bins away from the minimum solution and on the other hand, when a new item 
arrives, we want to repack as few items as possible. In the case of offline bin packing it is 
known that unless V = MV there is no polynomial time approximation algorithm for offline bin 
packing that produces a solution better than |OPr, where OPT is the minimum number of 
bins needed. For this reason, the most common way to deal with the inapproximability problem 
is the introduction of the asymptotic approximation ratio. The asymptotic approximation ratio 
for an algorithm A is defined to be limj;_i.oo sup{ Qpy^j) | OPT{I) = x}. This leads to the 
notion of asymptotic polynomial time approximation schemes (APTAS). Given an instance of 
size n and a parameter e E (0, 1], an APTAS has a running time of poly^n)-^^'^^ and asymptotic 



approximation ratio 1 + e, where / is an arbitrary function. An APTAS is called an asymptotic 
fully polynomial time approximation scheme (AFPTAS) if its running time is polynomial in 
n and K The first APTAS for offline bin packing was developed by Fernandez de la Vega & 
Lueker [13] , and Karmakar & Karp improved this result by giving an AFPTAS [19] (see survey 
on bin packing [B]). 

Since the introduction by Ullman of the classical online bin packing problem [23], there 
has been plenty of research (see survey [8]). The best known algorithm has an asymptotic 
competitive ratio of 1.58889 [21] compared to the optimum in the offline case, while the best 
known lower bound is 1.54037 [2]. Due to the relatively high lower bound of the classical online 
bin packing problem, there has been effort to extend the model with the purpose to obtain an 
improved competitive ratio. Gambosi et al. pH] presented a model where they allow repacking 
of items. They presented an algorithm which achieves ratio 1.33 by moving at most 7 items 
to a different bin each time a new item arrives. The idea of a dynamic packing was developed 
pretty early by Coffman, Garey and Johnson [4]. They developed and analyzed an algorithm 
for the dynamic bin packing model when arrival and departure of items are known in advance. 
Ivkovic and Lloyd [IT] presented an algorithm for dynamic bin packing having ratio 1.25. In 
their model items and their arrival and departure are not known in advance. The algorithm 
requires 0(logn) shifting moves, where a shifting move consists of moving a large item or a 
bundle of small items from one bin to another. In another work Ivkovic and Lloyd [16 gave an 
algorithm which achieves approximation ratio 1 + e by using amortized 0(logn) shifting moves. 
Concerning reassignment of jobs in scheduling, Albers and Hellwig [I] presented an algorithm 
for online scheduling with minimizing the makespan on m machines. The algorithm has a 
competitive ratio of am, where 02 = 4/3 and for m — )• 00 the competitive ratio Om converges 
to ~ 1.4659. The reassignment of jobs is bounded by 0{m). They also proved that it is not 
possible to obtain an algorithm with better competitive ratio than am without reassigning Q(n) 
jobs. 

The model we follow is the notion of robustness. Introduced by Sanders et al. [SU] it allows 
repacking of arbitrary items while the number of items that are being repacked is limited. To 
give a measure on how many items are allowed to be repacked Sanders et al. [20] defined the 
migration factor. It is defined by the complete size of all moved items divided by the size 
of the arriving one. An (A)PTAS is called robust if its migration factor is of the size f{j), 
where / is an arbitrary function that only depends on K Since the promising introduction of 
robustness, several robust algorithms have been developed. Sanders et al. [20] found a robust 
PTAS for the online scheduling problem on identical machines, where the goal is to minimize 
the makespan. The robust PTAS has constant but exponential migration factor 2*^*-e^°^ e^. 
In case of bin packing Epstein and Levin [TU] developed a robust APTAS for the classical bin 

packing problem with migration factor 2'^^^^'°^^^ and running time double exponential in K 
In addition they proved that there is no optimal online algorithm with a constant migration 
factor. Furthermore, Epstein and Levin [llj showed that the robust APTAS for bin packing can 
be generalized to packing d-dimensional cubes into a minimum number of unit cubes. Recently 
Epstein and Levin [T2] also designed a robust algorithm for preemptive online scheduling of jobs 
on identical machines, where the corresponding offline problem is polynomial solvable. They 
presented an algorithm with migration factor 1 — ^ that computes an optimal solution whenever 
a new item arrives. Skutella and Verschae |22j studied the problem of maximizing the minimum 
load given n jobs and m machines. They proved that there is no robust PTAS for this machine 
covering problem. On the positive side, they gave a robust PTAS for the machine covering 
problem in the case that migrations can be reserved for a later timestep. The algorithm has an 
amortized migration factor of 2'^'-e e^. 
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1.1 Our Results: 



An online algorithm is called fully robust if its migration factor is bounded by p(-), where p is a 
polynomial in K The purpose of this paper is to give methods to develop fully robust algorithms. 
In Section 2 we develop a theorem for a given linear program (LP) min{||x||^ \Ax > b,x > 0}. 
Given an approximate solution x' with value (1 + 6) LIN (where LIN is the minimum objective 
value of the LP) and a parameter a G (0, 5LIN], we prove the existence of an improved solution 
x" with value (1 + 6)LIN — a and distance \\x" — x'\\-^ < a{2/5 + 2). In addition, for a given 
fractional solution x' and corresponding integral solution y' , the existence of an improved integral 
solution y" with — = ( "+'" ) is shown (where m is the number of rows of A). Since 
both results are constructive, we propose also algorithms to compute such improved solutions. 
Previous robust online algorithms require an optimum solution of the corresponding ILP and 
use a sensitivity theorem by Cook et al. [7]. This results in an exponential migration factor in 
7 ( |m| [HI 1221 ) • In contrast to this we consider approximate solutions of the corresponding 
LP relaxations and are able to use the techniques above to improve the fractional and integral 
solutions. Furthermore we also prove an approximate version of a sensitivity theorem for LPs 
with modified right hand side b and b' . During the online algorithm the number of non-zero 
variables increases from step to step and would result in a large additive term. To avoid this 
we present algorithms in Section 3 to control the number of non-zero variables of the LP and 
ILP solutions. We can bound the number of non-zero variables and the additive term by 
0{eLIN) + 0{-^). In Section 4 we present the fully robust AFPTAS for the robust bin packing 
problem. We use a modified version of the clever rounding techniques of Epstein and Levin 
jlU| . This rounding technique is used to round the incoming items dynamically and control the 
number of item sizes. One difficulty is that we use approximate solutions of the LP. During the 
online algorithm items are rounded to different values and are shifted across different rounding 
groups. We show how to embed the rounded instance into another rounded instance that fulfills 
several invariants. By combining the dynanic rounding and the algorithm to get improved 
solutions of the LP and ILP, we are able to obtain a fully robust AFPTAS for the online bin 
packing problem. The algorithm has a migration factor of 0(l/e^) (or 0{\/e^) if the size of the 
arriving item is ^^(1)) and running time polynomial in ^ and t, where t is the number of arrived 
items. This resolves an open question of Epstein and Levin jlUj . We believe that our techniques 
can be used for other online problems like 2D strip packing, scheduling moldable tasks, resource 
constrained scheduling and multi-commodity flow problems to obtain online algorithms with 
low migration factors. 



2 Robustness of approximate LPs 

We consider a matrix A G W^^^ , a vector b G M^g ^-^d a cost vector c G M>g. The goal in 
a linear program (LP) is to find a x > with Ax > b such that the objective value c^x is 
minimal. We say x^^'^ is an optimal solution if c^x^^'^ = min |c"^x|^x > 6, x > o| and we 

define LIN = c^x^^"^ . In general we suppose that the objective function of a solution is positive 
and hence LIN > 0. We say x' is an approximate solution with approximation ratio (1 -|- 5) 
for some 6 G (0, 1] if \\x'\\i < (1 -|- 5)LIN . For the most part of the paper we will assume that 
C"^ = (1,1,..., 1) and therefore c^x^^"^ = x^^'^ ^ = LIN . The following theorem is central. 

Given an approximate solution x' , we want to improve its approximation by some constant. 
But to achieve robustness we have to maintain most parts of x' . We show that by changing x' 
by size of 0(f ), we can improve the approximation by a constant a. 

Theorem 1. Consider the LP Tai\i^c^x\Ax > 6, x > o| and an approximate solution x' with 
c^x' = (1 -|- 6)LIN for some 6 > 0. For every positive a < 6LIN there exists a solution x" 
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with objective value of at most c x" < (1 + 5)LIN — a and distance — < a{l/5 + 
1) " "^ Jv // = (1, 1, . . . , 1) then \\x' - < 2a(l/<5 + 1). 



Proof. We prove feasibility of the following LP 1 



Ax>b 
X > 

X > x' 



(LP 1) 



a{l/5 + l) 

T I 

C^ X 



c^x' 



c^x < {1 + 6) LIN 



a 



The assumption a < 6LIN implies c^x' = (1 + 6)LIN > a{l/6 + 1) and hence "^^t*^"!""*^^ < 1. 



Suppose that LP 1 is feasible and has a solution x" . Due to constraints 3 and 4 the distance 



between x" and x' can be bounded. Components with x[ > may be smaller compared to x' 



OPT 



by a(l/(5+l)^j^, while components x'/ with xf^'^ > may be larger than x^ by a(l/J+l) ^Vj./ • 



In the worst case, x' and x^^'^ have no common non-zero entries and hence, 



< 



«(i/<^ + i)(E.^li7 + Ei 



a(l/^ + l) " 



If c"^ = (1,1,...,!) then 



LP 1 



and - — T I 
feasibility of 
that each constraint of 
3 is fulfilled since x" - 
fulfilled since x" 



\\x'\\ 

< 1. Therefore a(l/5-|-l) ^ ^ , — < 2a{l/8 + l). It remains to prove 



We construct a solution x" by x" = (1- 



. We prove 



LP 1 



is satisfied for x". Note that x" > since "^^Jr^^^ < 1- Constraint 
X' — "^'^t'^"!""'''' x' + ^^^i^^^^^x*^^-^ > x' — "^^r'^"!""^^ x^ and constraint 4 is 
a;' _ ^^(i^x' + ^^^i^x^^"^ < x' + ^^^i^x^^^! Feasibihty for x" follows 

ai'i-/S+l) \ I 4^0PT "(1/^+1) 
T^/ ;-r^x , 



from Ax" = Ax'(l - ^^ii^) +^x 
objective value of x" is bounded by c'^x" 

1) + (Sgi^^^^ = c^^' - "(1/5 + 1) + 



cTa:'>a{l/5+l) 



T / 

C X 



> + = The 



T / 
C X 



a. 



□ 



From here on and for the rest of the paper we suppose that c^ = (l,l,...,l). 
Remark 1: 

Suppose x' has approximation ratio ||x'||^ = (1 + 5') LIN for some 6' > 0. By Theorem [T] the 
following LP is feasible if a' < 6' LIN. 

Ax>b 
X > 



x>x'-a{l/5' + l) 



\^ 111 

„OPT 



X' 



X < x' + a'(l/(5' + l) — 
{l + S')LIN-a 



Setting a' = Oijjjp^ for some 6 > yields feasibility for the following LP assuming \\x'\\^ 
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{1 + 6') LIN. 

Ax>b (LP *) 

X > 

x' 

x>x- aillb + V)-—— 
|p 111 

OPT 

X < x' + a(l/(5 + 1)^— 

If 111 

J2xi< {I + 5')LIN - a' 

Here, we use a' {1/5' + 1) = a^jp;^{l/5' + 1) = a{l/5 + 1). The condition that a' < 5' LIN is 
equivalent to the condition that ||2;'||j^ > a{\/5 + \) since ||a;'||;^ = {\ + 5')LIN and a{l/5 + l) = 
a'{l/5' + l). □ 
Remark 2: 

In many cases, we do not know the exact approximation ratio ||x'||-|^ = (1 + 5') LIN but the 
approximation guarantee \\x'\\-^ < (1 + 6) LIN for some 6 > 6' . Assuming \\x'\\i > a{l / 6 + 1) we 



can use feasibihty of LP * to prove the existence of a solution x" with ||x"||-^ < (1 + 5) LIN — a 



and ||x" — x'll;^ < 2a{l/6 + 1). The distance ||x" — x'||^ < 2a{l/5 + 1) follows again easily 

the ai 



LP 



from constraints 3 and 4 of LP * We derive the aimed objective value ||a;"||^ from the last 
constraint of 



||2;"||i < (1 + 6')LIN - a^f§j^ = (1 + 6)LIN - {6 - 5')LIN - a{l/5 + 

(l+5')L7iV>^(l/5+l) 

l)Tjhi ^ {l + 5)LIN-a{l/6 + l)^j^-a{l/S + l)j^ = {l + 6)LIN-a{l/6 + 

S'<5 

1) 1+5^ = (1 + 5) LIN — a j^pjf < (1 + 5) LIN — a. This proves, that it suffices to know an upper 
bound for the approximation to obtain an improved solution x". □ 
Of course, one major application of Theorem[T]is to improve the approximation. But we can 
also apply Theorem[l]to obtain a variant of the theorem of Cook et al. [7J for the sensitivity anal- 
ysis of an LP. Consider the following problem: Let x' be a solution of min { ||x|| ^ \Ax > b' ,x > 0}. 
Find a solution x" for LIN 2 = min {||2;||-|^ \Ax >b",x>0} with changed right hand side such 
that ||x" — x'll]^ is small. A theorem of Cook et al. [7] states that there exists a x" satisfying the 
LP and H^;" — ic'||^ < nA \\b" — b'\\^, where A is the largest subdeterminant of A. This result 
is not satisfying if A and n are too big, especially if they are exponential in m. By letting loose 
of optimal solutions we obtain a corollary that is much more appropriate to derive fully robust 
algorithms. In contrary to the theorem of Cook et al. [7] the amount of change in the solution 
does not depend on the determinant nor on the dimensions of A but on the approximation ratio 
of the solution. 

Corollary 2. Consider the linear program LP^ defined by min{||x||j^ \Ax >b',x>0} and an 
approximate solution x' with \\x'\\^ < {1 + 6)LINi (6 > 0) and \\x'\\-^ > a{l/5 + 1). There exists 
a solution x" of LP^ defined by mm{\\x\\^ \Ax > b" , x > 0} with ||x"||^ < (1 + 6) LIN 2 such that 



'2 

^ b"-b'- 
components ' ' . 



the distance ||x" — x'||^ < (| + 7) ^ ^ ^ ^ where Ci = maxj Aij and " ^ " is a vector having 



Proof. Suppose there is only one index i where b'^ / 6". Consider the 2 cases: 

b" —b' 

Case 1: 5^ < b'J We increase x'j by ' , where j is the index with the maximum entry in row i. 
This way we make sure that the so modified x' covers the larger b'- since now {Ax')i = b'-. Since 
we simply increase x' to cover the larger b" we may worsen the approximation by an additive 
term of at most \, ' . 

Case 2: b'^ > b'J. In this case we do not modify component i of x', but since a smaller b'- has 
to be covered the optimal value of a solution may decrease. Let LINi be the optimal value of 
LP^ and let LIN2 be the optimal value of LP^. The inequahty LIN2 < LINi leads to 
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a contradiction since we can increase an optimal solution x^^^^ of LP^ to cover the lager h' like 
we did in case 1. Modifying x^^^^ this way would lead to a smaller optimal solution. Therefore 
the optimal solution of LP can not decline by more than ' ' . Using approximate 



solution for LP^ yields therefore \\x'\\^<{l + S){LIN2 + = (1 + 5) LIN 2 + (1 + 5)-^. 

Iterating over all components 1 < i < m and changing the solution according to the cases 



would result in an approximate solution of at most {\ + 5)LIN + {1 + 5) 



. Using Remark 



b"~b' 



guarantees the existence of a solution x for LP^ having value 



b"-b' 



2 with a = (1 + (5) 

(1 + 6)LIN - a and ||x - x'\\^ < {2/6 + 2)a = {2/6 + 2)(1 + 6) 
Modifying x according to the cases yields therefore a solution x" with 
the distance between j;" and we get \\x" ^ --'n ' ^ 



X \\i < \\x 



< {2/6 + 6) 
x''\\i < {l + 6)LIN2. For 



b"-b' 



< 



+ 7) 



□ 



Note that if A is an integral matrix without zero rows, each component q is at least 1. 



3 Algorithmic Use 

Let x' be an approximate solution of the LP with < (1 + 6)LIN. In Theorem [Tj we have 

proven the existence of a solution x" near x' with ||a;"||;^ < (1 + 6) LIN — a. We are looking 
now for algorithmic ways to calculate this improved solution x" . We present two algorithms 



that basically rely on solving an LP. According to LP * we sp lit x' i nto a fixed part x-^*^' and 
a variable part x'"'^'^. The variable part is defined according to LP * by x^"'' = ^^^jj^Tjp^^;' and 

the fixed part by x-^*^ = x' — x^"""^. By assigning the variable part in a better way, Theorem [l] 
and Remark 2 state under the assumption that > a{l / 6 + 1) that we can improve the 

objective value by a. We denote with 6^°'' = b — A{x-I'^^) the part which has to be reassigned. 
The algorithm works as follows: 

Algorithm 1. 

1. Set x'"'"'' := ^'jl^x', xf'"" := x' - x™*" and b""'' := b - A(x-^^^) 

2. Solve the LP x = min{||x||i \Ax > x > 0} 

3. Generate a new solution x" = x-^*^ + x 

If X is a basic feasible solution, compared to x', our new solution x" has up to m additional 
non-zero components. 

Theorem 3. Given solution x' with \\x'\\-^ < (1 + 6)LIN and ||x'||^ > a{l/6 + 1). Algorithm^ 
returns a feasible solution x" with ||x"||^ < (1 + 6)LIN — a and the distance between x' and x" 
is \\x" - x'W^ < 2a{l/6 + l). 

Proof. Solution x" is feasible because A{x") = A(x^*^+x) = ^(x-''^^)+^(x) > A{xf^'')+b'"''^ = b. 

For the approximation we use Remark 2, which guarantees the existence of a solution with 

objective value < (1 + 6) LIN — a by leaving the part x-^^^ o f the s olution x' unchanged. The 

unchanged part X'^*^ is defined by using the lower bounds of LP *, x" > x' — a{l/6 + 1) Zu ■ 

I I 11^' 111 

Placing X optimally leads therefore to the aimed approximation. Since the ||x^'^''||;^ and ||x||^ 
are bounded by a{l/6 + 1) the worst possible distance between x' and x" is ||x" — x'||j^ = 

|(x-''^^ + x)-(x^^^ + x^'^0||^ = ||x-x^'^''||i < Pill + llx^'^nii <2a(l/5 + l). □ 

In Algorithm [T] we use an optimal LP solver as a subroutine. In many cases, like for 
example bin packing, the corresponding LP relaxation is hard to solve and the running time 
for computing an optimal solution is very high. For the following algorithm it is sufficient to 
compute the LP approximately, which in general can be performed more efficiently. We assume 
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that \\x'\\^ > 2a{l/5 + 1) because the double amount has to be reassigned to achieve the same 
improvement in the approximation as in Algorithm [T] 



Algorithm 2. 

1. Set x'"^'' :-- 

2. Solve X = 



2a(l/<5+l) / fi^ ._ 



mm 



kill \Ax > 6^"'' 



3. If 



fix 



+ X 



< x' 



set X ' 



x' - x""""- and b'""'' := b - 

X > 0} approximately with ratio (1 + 6/2) 



Theorem 4. Let x' be a solution with \\x'\\-^ < (1 + S)LIN and \\x'\\-^ > 2a{l/6 + 1). Then 
Algorithm^ returns a feasible solution x" with approximation guarantee (1 + 5) LIN — a and 
||x"-x'||i <4a{l/S + l). 

Proof. The property that ||rE" — x'H-^ < 4a(l/(5 + l) follows by Theorem 3 and the fact, that x"""""^ 
has the double size 2a{l/S + 1) compared to x""""^ defined in Algorithm [T| Furthermore we have 
to show that at the end of the algorithm ||x"||]^ < (1 + S)LIN — a. Suppose ||x'|| = (1 + 6')LIN 
for some 5' < 6. Using the assumption 2a{l/6 + 1) < H^'H^ < (1 + 6) LIN implies that 
2a < ^^ly'^^^^j^ = 5LIN . Consider the case that 25' < 5. In this case x" has the aimed 
approximation since \\x"\\^ < \\x'\\^ = (1 + 8')LIN < (1 + 5)LIN - 5/2LIN < (1 + d)LIN - a 
using 2a < 5LIN . Thus in the following we assume 5 < 26' . Suppose we solve the LP in step 



2 optimally, 
feasibility of 



n this case, Algorithm [2] is identical to Algorithm [T] using improvement of 2a. By 



LP 



we know there exists a solution x" with ||x"|| < (1 + 6')LIN — 2a{ 



1/5+1 
l/S'+l' 



x 



OPT 



This implies, that an optimal solution x^^'^ of the LP min{||x||]^ \Ax > b'"°'^,x > 0} is of size 

2a{l/S + 1) — 2a Ijgr^i ■ Solving the LP approximately with 

ratio (1 + 6/2), solution x has an additional term 5/2 



< ||x^ 



2a V-^+i 



%OPT 



bounded by ||x||;^ < 



,OPT 



+ 6/2 



x 



OPT 



< W 



this results in the approximation for x-^*^ + x as follows. 



The value of \\x\\-^ is therefore 
-2a^ + a(l + 5)-a(l^). Finally 



x^'^ + X 



Ix^'^HL + ll^l 



=(1 + 6')LIN - 2a(^^|^) + a{l + 6)-a{-^^^ 



■1/6' + 1' 



■1/6' + 1' 



--{1 + 6)LIN -{6- 6')LIN - 2a(^/|±i) + a{l + 6)- a(-^-^) 



1/6' + 1' 



LIN'>2a/5 § 

< {1 + 6)LIN -2a{—^) + a{l + 6)-2a{ 



6'. 



[1 + 5)LIN + a{- 



-25 + 26' + 6 + 5'^. 



1/6 + 1 
1/6 + 1. 



■1/6' + 1 
1 + 6 



a[ 



=(1 + 5)LIN 
--{1 + 6)LIN 
--{1 + 6)LIN 



a + a{ 
a + a{ 
a + a{ 



6 

26' + 6"^, 



2a{ 



a{ 



1/5' + 1 
1 + 6 



2a{ 



1/6' + 1' '1/6' + 1 
1/(5 + 1, , 1 + 6 



) 



6 ' '1/6' + 1' 
2 + 26' + 5'^/6' + - 2 



25 -(5 -(52, 



(5(1/(5' + 1) 
2(5' + (52/<5'-35, 



6{l/6> + 1) 



-) 



=(1 + 5)LIN - a + a 
<(1 + 5)LIN 



X5-5'){-2 + l,) 
(5(1/(5' + 1) ' 



a 
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The last inequality holds because a{ ^^ ^(i/^'+T/' ) 
By the last step of the algorithm we know that 
5)LIN-a. 



< since 6-6' > and -2+^ < <^ 6 < 26'. 



< 



X 



fix 



+ X 



and thus 



kill < (1 + 

□ 



In some cases we may not want to get a guaranteed approximation, but a guarantee that 
our solution x' is getting smaller by some a. This works if the approximation ratio of x' is worse 
than (1 + 6). The following corollary states, that if we use Algorithm |2] on a solution x' with 
= (1 + 6') LIN for some 6' > 6 the objective function of our new solution x" decreases by 
at least a. 



Corollary 5. Let \\x'\\-^ = {1 + 6')LIN for some 6' > 6 and \\x'\\^ > 2a{\/6+\). Then Algorithm 
returns a solution x" with < — a = (1 + 6') LIN — a and \\x" — x'\\^ < 4a(l/(^ + 1). 

Proof. Suppose like in the proof of Theorem 4 that we solve the LP in step 2 optimally. In this 
ca se, A' 
of 



LP 



.gorithm|2]is identical to Algorithm[T]using improvement of 2a and therefore by feasibility 
e know that it returns a solution x" with ||x"|| < {1 + 6')LIN — 2a{Yjjrp[). An optimal 



-,OPT 



< W 



solution x*^^"^ of the LP min{||x||^ \Ax > b"°'^\x > 0} is therefore of size 
2a Ij^r^i = 2a(l/6 + 1) — 2a Ij^r^i ■ Since we actually solve the LP approximately with ratio 
{1 + 6/2), solution x has an additional term of 5/2 

by \\x\\i < - '^a^^l±^ ^ - -^^+^ 



,OPT 



and the value is therefore bounded 



^ — 2ajjy:p^ + a{l + 6) — a{ ij^i+i ) according to the proof of Theorem 5. By 



construction of x" we get < 



fix 



+ X 



< {1 + 6')LIN -a. 



Since 6' > 6 we know that —2a{ Ij^i^i ) ^ ~2q; and that a(l + 6) — < Hence 



[l + 6')LIN-2a{^^) + a{l + 6)-a{jj±^^). 



□ 



3.1 Integer Programming 

In this section we discuss how we can apply results from the previous sections to integer pro- 
gramming. Consider a fractional solution x' of the LP and a corresponding integral solution 
y' . By rounding each component x'^ up to the next integer value, it is easy to get a feasible 
integer solution y' with an additional additive term \\y'\\i < + C, where C is the number 

of non-zero components. We can apply any of the previous algorithms to x' to get an improved 
solution x" . But our actual goal is to find a corresponding integer solution y" with improved 
objective value < (1 -|- 6) LIN -\- C — a such that the distance between y" and y' is small. 

In the following we present two algorithms that compute a suitable y" with improved objective 
value and small distance between y" and y'. Note that the straight forward approach to simply 
round up each component x'- leads to a distance between y" and y' that depends on C and hence 
(depending on the LP) is too high. Designing the algorithms, there seems to be some trade 
off between the number of non-zero components and the distance between the integer solutions 
y' and y". The first algorithms tries to minimize the distance between y' and y" while the 
second guarantees better approximation of \\y'\\i and \\y"\\i while the distance between them 
increases. The existence of an algorithm combining both good properties, low distance and 
good approximation guarantee of y' and y", is an interesting question. 

In Algorithm [s] we focus on how much components of x' need to be reduced to achieve the 
improved approximation guarantee. This defines the migration factor in robust bin packing. 
The actual worst case distance between y" and y' is larger and however can only be bounded by 
0{m + 1/6). Like in the previous algorithms, we assume that > a{l/6 + 1). We require 

x' to be a solution with approximation guarantee < (1 -|- 6) LIN and we require y' to be 

an integer solution with approximation guarantee \\y'\\i < (1 -|- 6) LIN + n. For every 1 <i <n 
we suppose that x'^ < y'^. For a vector z G K>0) let ^(-2) be the set of all integral vectors 
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V = {vi, . . . ,Vn)'^ such that < Vi < zi. Given LP solution x' and integer solution y' with the 
described properties above. The algorithm performs in the following way. 

Algorithm 3. 

1. If possible choose vector c G V{y' — x') with ||c||^ = a and return y" = y' — c and x" = x' . 
Otherwise choose c G V{y' — x') such that \\c\\-^ < a is maximal. 

Set y = y' — c. 

2. Set x'"^'' := "^jg+^^ rcS x^^^ := x' - x""'' and h""'' := b - ^(x^^^) 

3. Compute an optimal solution x of the LP min{||x||j^ \Ax > b'"^'^,x > 0} 

4. Set x" = x^^'-" + X 

5. For each 1 < i < n set iji = max{ [x"] , yi} 

6. If possible choose d G V(jj—x") such that \\d\\i = a(l/5+l) otherwise choose d G V{y—x") 
such that \\d\\^ < a{l/5 + 1) is maximal. 

7. Return y" = y — d. 

Theorem 6. Let x' be a solution of the LP with \\x'\\-^ < (1 + 5)LIN and \\x'\\^ > a{l/5 + 1). 
Let y' be an integral solution of the LP with \\y'\\i < (1 + 6)LIN + n where y[ > x[ for each 
i = 1, . . . ,n. Then Algorithm^ returns an integral solution y" with ||y"||;^ < (1 + 5)LIN + n — a 
such that J2y'.' Ky'Xy'i - y'i) < a{l/S + 2). 

Proof. Feasibility: Feasibility for x" and approximation ||x"||]^ < (1 + 5) LIN — a follows from 
Theorem 3. Step 2,3 and 4 are identical to Algorithm [T} Feasibility for the integer solution 
y" follows from the fact, that for every component i we have y'/ = jji — di > x'l and hence 
Ay" > Ax" > b. 

Size of reduction of y': The only steps where components of y' are changed are in step 1, 5 
and 6. In step 1 we change y' to obtain y, in step 5 we change y to obtain y and in step 6 we 
change y to obtain y". Summing up the change in each step leads therefore to the maximum 
possible size of reduction of y' compared to y" . In step 1 there are c < a components of y' 
which are being reduced. In step 5 no components of y are being reduced and in step 6 there 
are d < a{l/6 + 1) components of y which are being reduced to obtain y" . Hence there are at 
most a{l/5 + 2) components of y' which are being reduced to obtain y" . 
Approximation: It remains to prove, that y" has approximation ratio (1 + 5)LIN + n — a. 
Case 1, ||c||-|^ = a: In this case, the algorithm returns in step 1 solution y" = y with = 
\\y'\\i — a and the algorithm terminates. Otherwise, if ||c||;^ < a we have for every component i, 
that yj — x^ < 1 and ||x"||-^ < (1 + 6)LIN — a. Note that steps 2-4 are equivalent to Algorithm 

m 

Case 2, \\d\\-^ < a{l/6 + 1): In this case y" — x" = iji — di — x" < 1 for i = 1, . . . , n, since \\d\\i 
is chosen maximally. Using ||x"||-|^ < (1 + 5) LIN — a and y" < x" + 1 for i = 1, . . . , n we have 
||y"|li < {1 + S)LIN + n-a. 

Case 3, = a{l/6 + 1): Let rh be the number of components with x" > yi. Next we 

compare the vector y with x". Using x" > x-^*^ and the definition of y in step 5 we obtain 
\\y - 111 = Y.x'^<yM - O + Ex;'>j;,(r<l - O < Y.x'^<yM - 4 ) + ™- The fact that 
fji — x'^ < 1 for i = 1, . . . , n and ||x'||-|^ — x-^*^' ^ = ||x''"'"|| < a(l/(5 + 1) and the fact there are at 

most n — rh components with x" < yi yield that J2x"<yi{yi ~ ^f") = J2x"<yi{yi ~ + ^i"'^) ^ 
n — rh + J2x"<y' ^r"^ < n — rh + a{l/6 + 1). As a result we can bound ||y — x"\\-^ < J2x"<y'iyi ~ 
x{^^) + fn < n + a{l/5 + 1). Since y" = y — d and = a{l/6 + 1), our integer solution y" has 
the aimed approximation guarantee of \\y"\\i = \\y\\i — \\d\\i < ||x"|| + a{l/5 + 1) + n — \\d\\i = 
\\x"\\+n<{l + 5)LIN + n-a. □ 
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The running time of the above algorithm depends on the number of non-zero components 
and the time to compute an optimal solution of an LP. The algorithms computes an integral 
solution y" with < (1 + 5)LIN + n — a for given fractional and integral solution. In 

many cases, like bin packing, the dimension n is very large and provides thus a large additive 
term in the approximation. The following algorithm describes how this large additive term can 
be avoided. On the other hand the difference between y' and y" increases to Let 
x' be an approximate solution of the LP min{||x||]^ \Ax > b, x > 0} with \\x'\\i < (1 + 6)LIN 
and > a{l/5 + 1). Furthermore let y' be an approximate integer solution of the LP with 

lly'lli < (1 + 25)LIN and \\y'\\^ >{m + l){l/5 + 2) and y[ > x\ for i = 1, . . . , n. In addition we 
suppose that both x' and y' have exactly K < 6 LIN non-zero components. Our goal is now to 
compute a fractional solution x" and and integer solution y" having improved approximation 
properties and still < SLIN non-zero components. For a vector z £ K>0) let ^{^) be the set 
of all integral vectors v = (ui, . . . , Vn)^ such that < Vi < Zi. Furthermore we denote with 
oi , . . . , ax the indices of the non-zero components y'^. such that y^^ < • • • < y'a^^ are sorted in 
non-decreasing order. 

Algorithm 4. 

1. Choose I maximally such that the sum of smallest i components !,...,£ is J2i<i<iy'ai ^ 
(m + l)(l/(5 + 2) 



x'i if i = aj for j < e f if i = aj for j < 




2. Set xr = { aii/5+i) ^, ' andm 

3. Set X'^*^ = x' — x""'^ , b""""^ = b — A{x-^'^^) and compute an optimal solution x of the LP 
min{||x||-^ \Ax > > 0} 

4. Set x" = x^*^ + X 

5. For each 1 < i < n set jji = max{[x"] , yi} 

6. If possible choose d G V(jj — x") such that \\d\\i = a{l/5+l) otherwise choose d G V{jj—x") 
such that \\d\\^ < a{l/6+ 1) is maximal. 

7. Return y" = y — d 

Theorem 7. Let x' be a solution of the LP with \\x'\\i < (1 + 5)LIN and \\x'\\^ > a{l/S + 1). 
Let y' be an integral solution of the LP with \\y'\\i < {1 + 25)LIN and \\y'\\i > {m + l)(l/5 + 2). 
Solutions x' and y' have both exactly K non-zero components and for each component we have 
x'l < y[- Then Algorithm^ returns a fractional solution x" with < (1 + 6)LIN — a and 

an integral solution y" with \\y"\\i < (1 + 25) LIN — a. Both x" and y" have the same number 
of non-zero components with x'- < y'- and the number of non-zero components is bounded by 
6 LIN. The distance between y" and y' is bounded by ||y" — = 

Proof. Feasibility: Feasibility and approximation for the fractional solution x" follow easily 
from correctness of Algorithm [T] and the fact that removing additional components x'^^ , • • • , x'^^ 
and reassigning them optimally does not worsen the approximation. Each integral component 
yi is by definition (step 5) greater or equal than x'-. By choice of d step 6 and 7 retain this 
property for y" and imply thus feasibility for y" . 

Distance between y" and y': The only steps where components of y' are changed are step 
2, 5 and 7. In step 2 we change y' to obtain y, in step 5 we change y to obtain y and in step 7 
we change y to obtain y". Summing up the change in each step leads therefore to the maximum 
possible distance between y" and y' . In step 2 of the algorithm £ components of y' are set to 
zero to obtain y, which by the definition of i results in a change of at most (m -\- 1) {1 / 5 -\- 2) . 
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We define L by L = J2i<i<£ v'a^ with < L < (m + l)(l/5 + 2). In step 5, tlie only components 
yi being changed are the ones where x'l is larger than yj. So the change in step 5 is bounded 
by Y.x'|>n'^\x'l^ - Vi) = Y.x'l>n^\4''' + - m) < E<>j;,(ra;f''l -yi + [xi]) < Y.x'l>y,\^i\ 
by knowing that — y « < since = yi = ^ \i i = aj for a j < £ or \x{^^^ < [x^] < y-. 

Furthermore we can bound J2x">yi\^i\ — ll^lli + < ll^^^^lli + since x is a basic feasible 
solution and can be bounded by L + a{l/6 + 1) (i.e. we get L for the size of components 



X' 



ai) • • ■ '^ciK P^^^ ^i>t "^IMI^^^ — a(l/(5 + 1) for the remaining ones). Therefore we have 
||y — y|li < i + a(l/(^ + 1) + m. In step 7, ||y" — y^^ — Mlli ^ a(l/(5 + 2) + m. In sum this 
makes a total change of at most (m + + 2) + L + aiXjb + 1) + m + a(l/(5 + 2) + m < 

2(m + l)(l/(5 + 2) + 2m + a(2/(5 + 3) = OiJ^). 

Number of components: The property that x' and y' have the same number of non-zero 
components together with the property that y\ > x'^ implies that x'^ > whenever y| > 0. This 
property holds also for x-^^^ and y since a component yi is set to zero if and only if xp^ = 0. 
Notice that y" = y — d > x" . Suppose by contradiction that there is a component i with x" = 
and > 0, then jji = y'- + di > and by definition of y we obtain yi > 0. In this case we have 
x{^^ > 0, which gives a contradiction to x" = = x{^^ + Xi > 0. Using the property that x" 
and y" have the same number of non-zero components, it is sufficient to prove that the number 
of non-zero components of x" is limited by 6LIN. Our new solution x" is composed of x-^"' and 
X. Solution x^^^ has K — £ non-zero components, since in step 2 we set £ components of x^^^ 
to zero. Being a basic feasible solution, x has at most m non-zero components and hence x" 
has at most K + m — £ non-zero components, li £ > m, then x" has < K < 6LIN non-zero 
components. So let i < m: The total number of non-zero components after step 4 is {K +m — £). 
We now prove that this number is bounded by 5LIN. Parameter £ is chosen to be maximal, 
therefore J2i<i+i v'ai ^ [m + l){l/5 + 2). Hence, the average size of components y'^_^, . . . ^y'a^j^-^ 

is greater than ("^+^^)07'^+^) > ^""'"^^^^^^^^ > 1/^ + 2. Since the components are sorted in 
non-decreasing order, every component y[ with i > i + 1 has size > 1/5-1-2. Summing over 
all non-zero components of y' yields the following inequality: \\y'\\i = J2iLe+2 Uai + Vai+i + ^ ^ 
{K-£-l){l/6 + 2) + y'^^^^+L > {K -£-l){l/6 + 2) + {m + l){l/6 + 2) = {K -£+m){l/6 + 2). 
Using that \\y'\\^ < (1 -h2(^)L/7V yields {1 + 25)LIN > {K - £ + m){l/S + 2). Dividing both sides 
by (1/6 + 2) gives {K — £ + m) < 5LIN. This shows that the number of non-zero components 
of x" and y" is at most 6LIN. 

Approximation: Casel: = a(l/(5 -|- 2) -|- m 
The following inequalities ||y||-^ < \\y\\i + L + a{l/5 + 2) + m = \\y'\\i + a{l/5 + 1) + m and 
lly'lli ^ (1 + 25)LIN together yield the aimed approximation ||y"||i = \\y\\i — \\d\\i = \\y\\i — 
a{l/6 + 2) - m < (1 + 26)LIN - a. 
Case2: \\d\\-^ < a{l/6 + 2) + m 

Since d is chosen maximally, y'/ — x" < 1 for every components i = 1, . . . ,n. Since ||x"||-|^ < 
(1 + 5) LIN — a and y" has at most 5 LIN non-zero components ||y"||i is bounded by {1 + 5)LIN — 
a + 6LIN = {1 + 25)LIN -a. □ 

Instead of using an optimal LP solution in Algorithm [3] and |4] we can solve the LP approxi- 
mately with a ratio of (1-1-5/2). The following algorithm is basically a combination of Algorithm 
[2] and Algorithm |4j We could also combine Algorithm [2] and Algorithm |3] to obtain similar re- 
sults. We make the following assumption for the fractional solution x' and the corresponding 
integer solution y': Let x' be an approximate solution of the LP min {||x||]^ \Ax > b,x > 0} with 
< (1 + 5)LIN and \\x'\\-^ > 2a{l/5 + 1). Let y' be an approximate integer solution of 
the LP with \\y'\\-^ < LIN + 2C for some value C > 6LIN and with \\y'\\^ > {m + 2){l/6 + 2). 
Suppose that both x' and y' have only K < C non-zero components. For every component i 
we suppose that y'^ > x[. Furthermore we are given indices ai, . . . , ax, such that the non-zero 
components y'^. are sorted in non-decreasing order i.e. y'^_^ < . . . < VaK- 
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Algorithm 5. 

1. Set x""'' := 2^ii^^x', x^"' := x' - x""""^ and 6™'' = 6 - A(x-^") 

2. Compute an approximate solution x of the LP min{||j;||^ \ Ax > 6^"',x > 0} with ratio 
(1 + 5/2) 

3. If x-^*^ + X ^ > \\x'\\-^ then set x" = x' , y = y' and goto step 9 

4- Choose the largest £ such that the sum of smallest components y'l, ■ ■ ■ ^y[ is J2i<i<eya — 
{m + 2){l/S + 2) 

-fix I if i = o-j for j < £ _ I if i = ttj for j < £ 

5. For all i set xf = < r-^. and yi = \ 

\xl else yy[ else 

6. Set x = X + x'^ where x'^ is a vector consisting of components Xa^ , • • • , ■ Reduce the 
number of non-zero components to at most m + 1. 

7. x" = x-^*^ + X 

8. For all non-zero components i set yi = max{ [x"] , yj} 

9. If possible choose d £ V{y—x") such that \\d\\-^ = a{l/6+l) otherwise choose d G V{y—x") 
such that \\d\\^ < a{l/5 + 1) is maximal. 

10. Return y" = y — d 

Step 6 of the algorithm can be performed using a standard technique presented for example in 
[3]. Arbitrary many components of x can be reduced to m+1 without making the approximation 
guarantee worse. We formulate the following theorem and corollary such that we can directly 
use it in the next section. 

Theorems. Letx' be a solution of the LP with \\x'\\^ < {1+6)LIN and\\x'\\^ > 2a(l/5+l). Let 
y' be an integral solution of the LP with \\y'\\i < LIN + 2C for some value C > 6 LIN and with 
Wv'Wi ^ (m + 2)(l/5 + 2). Solutions x' and y' have the same number of non-zero components and 
for each component we have x[ <y[. The number of non-zero components of x' and y' is K with 
K <C . Then Algorithm^ returns a fractional solution x" with ||x"||;^ < (1 + S)LIN — a and an 
integral solution y" where one of the two properties hold: \\y"\\i = a or = ||x"||-^+C. 

Both, x" and y" have at most C non-zero components and the distance between y" and y' is 
bounded by \\y" - y'\\^ = 0{^). 

Proof. Note that the first 3 steps are equivalent to Algorithm [2| In steps 4-6 the number of 
non-zero components x^^ , . . • , x^^ are reduced. As we apply a method that does not increase 
the objective value we obtain by Theorem 4 that ||x"||-^ < (1 -|- 5)LIN — a. Steps 4-9 are similar 
to Algorithm 4| The main difference is that components x^^ , • • . , x^^ are not assigned by the 
LP but are added to the LP solution afterwards in step 7. 

Distance between y" and y': As in Theorem 7, the steps where components of y' are changed 
are steps 5,8 and 10. By definition of £ the change of y' in step 5 is bounded by (m + 2)(l/6 + 2). 
As shown, the change in step 8 is bounded by J2x">yi T^f*^ + ^^l ~ Vi) — J2x">yi — II 111 
and ||[xj]||]^ < 2a(l/(5 -|- 1) -|- L -|- 1, whereL = J2i<i<eyai- The change in step 10 is bounded 
by IMIli ^ 2a(l/(5 + 2) + m + 1. Therefore the total change between y' and y" is bounded by 

Number of components: According to Theorem 7, the number of nonzero components of y" 
is equal to the number of non-zero components of x" which equals K — £ + m-\-l (the number of 
non-zero components of x is bounded by m + 1). We distinguish between the two cases where 
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I > m + 1 and I < m + 1. In the case where £ > m + 1 the number of components of x" is 
smaller than K and hence bounded by C. Consider the case where i < m + 1. By definition of 
i we know that L + y'^_^_^ > (m + 2)(l/5 + 2). Using the argument in the proof of Theorem 7, we 
obtain the following inequality: \\y'\\i = Y,i=e+2 Ui + y't+i+L = {K — i—l){l/5 + 2)+y'^^-^+L > 
{K - £ - l){l/6 + 2) + {m + 2){l/6 + 2) = {K -£ + m + l){l/6 + 2) Using that \\y'\\^ < LIN + 2C 

yields LIN + 2C>{K-i + m + l){l/5 + 2). As = ^-^f^^ ggf = C we 

obtain that {K - £ + m + 1) < \i^s^^ ^ ^■ 

Approximation: According to Theorem 7 we distinguish between the two cases where \\d\\-^ = 
2a{l/6 + 2) + m + 1 and < 2a{l/6 + 2) + m + 1. In the second case where < 

2a{l/6 + 2) + m + 1 we know that is bounded by plus the number of non-zero 

components of x" since d is chosen maximally. Hence ||y"||x < + C. In the case where 

\\d\\^ = 2a{l/6 + 2) + m + 1, we know \\y"\\-^ < \\y\\-^ - 2a{l/6 + 2) - m - 1. As \\y\\-^ < 
\\y\\i + II 111 ^ lly'lli + 2a{l/5 + 1) + m + 1 we get ||y"||;^ < ||y'||i — a- Note that we can also 
make the general claim for y" that ||y"||x < LIN + 2C — a. □ 

The following corollary is an analog to Corollary [5] which states what Algorithm [5] is doing 
if the approximation ratio of x' is worse than (1 + 6). We will need this corollary in the next 
section as we have no true control about the approximation ratio of x' . During the bin packing 
algorithm new columns might appear in the LP, which might change the optimal solution and 
therefore the approximation ratio of a solution x' . 

Corollary 9. Let \\x'\\-^ = (1 + S')LIN for some 6' > 5 and \\x'\\^ > 2a{l/6 + 1) and let 
lly'lli < LIN + 2C for some C > 6' LIN and \\y'\\^ > (m + 2){l/6 + 2). Solutions x' and y' 
have the same number of non-zero components and for each component we have x[ < y[. The 
number of non-zero components of x' and y' is K with K < C . Then Algorithm returns a 
fractional solution x" with ||x"||-|^ < ||x'||-|^ — a = (1 + 5')LIN — a and an integral solution y" 
where one of the two properties holds: \\y"\\i = ||y'||i — a or \\y"\\i = \\x'\\i — a + C. Both x" 
and y" have at most C non-zero components and the distance between y" and y' is bounded by 
\\y" -y'\\^=0{^^). 

Proof. Note that steps 1-3 are basically identical to Algorithm [2| Hence Algorithm [5] returns by 
Corollary 6 a fractional solution x" with ||2;"||-|^ < ||x'||;^ — a. The distance between the integral 
solutions y' and y" are independent of the approximation ratio of x' . Hence the distance between 
y' and y" is according to Theorem [s] bounded by 0{ ^^°' ). The number of non-zero components 
of x" and y" is by the proof of Theoremjsjbounded by the number K < C oi non-zero components 
of y' or by ^y^^2 — '^^^ approximation guarantee for y", that ||y"||]^ < ||?/'||i — ol follows if 
= 2a{l/5+2)+m+l. If \\d\\^ < 2a{l/5+2)+m+l then \\y"\\^ < \\x"\\^+C < \\x'\\^+C-a. 
We can also make the general claim for y" that ||y"||^ < ||y'||]^ — a. □ 



4 AFPTAS for robust bin packing 

The goal of this section is to give a fully robust AFPTAS for the bin packing problem using the 
methods developed in the previous section. For that purpose we show at first the common way 
how one can formulate a rounded instance of bin packing as an ILP. In Section 4.2 we present 
abstract properties of a rounding that need to be fulfilled to obtain a suitable rounding and in 
Section 4.3 we present the used dynamic rounding algorithm. The crucial part however is the 
analysis of the dynamic rounding in combination with ILP techniques. Since the ILP and its 
optimal value are in constant change due to the dynamic rounding, it is difficult to to give a 
bound for the approximation. Based on the abstract properties we therefore develop techniques 
how to view and analyze the problem clS Si whole. 
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The online bin packing problem is defined as follows: Let It = {ii, ■ ■ - it} be an instance with 
t items at time step t G N and let s : — ^ (0, 1] be a mapping that defines the sizes of the 
items. Our objective is to find a function Bt : {ii, ■ ■ ■ ,it} — >■ N"*", such that J2i:Btii)=j — ^ 
for all j and minimal maxi {Bt{i)} (i.e. Bt describes a packing of the items into a minimum 
number of bins). We allow to move few items when creating a new solution -Bt+i for instance 
/f+i = ItD{it+i}. Sanders et al. ^0] and also Epstein and Levin [TU] defined the migration factor 
to give a measure for the amount of repacking. The migration factor is defined as the total size of 
all items that are moved between the solutions divided by the size of the arriving item. Formally 
the migration factor of two packings Bt and Bt+i is defined by J2j<t:Bt(ij)j^Bt+i{i.j) ^ih)/ ^('''t+i)- 

4.1 LP-Formulation 

Let I be an instance of bin packing with m different item sizes si, . . . , Sm- Suppose that for 
each item £ I there is a size sj with s(ifc) = Sj. A configuration Cj is a multiset of sizes 
{a{Ci, 1) : si,a{Ci,2) : S2, • • .a{Ci,m) : Sm} with Ei<j<m < where a{Ci,j) denotes 

how often size sj appears in configuration Cj. We denote by C the set of all configurations. Let 
\C\ = n. We consider the following LP relaxation of the bin packing problem: 

min \\x\\-^ 

Xia{Ci,j) > bj VI < j < m 

Xi>0 VI < j < n 

Component bj states the number of items z in I with s{i) = Sj for j = 1, . . . ,m. This LP- 
formulation was first described by Eisemann [9J. Suppose that each size sj is larger or equal to 
e/2 for some e G (0, 1/2]. Since the number of different item sizes is m, the number of feasible 
packings for a bin is bounded by \C\ = n < (| + 1)™. Obviously an optimal integral solution 
of the LP gives a solution to our bin packing problem. We denote by OPT(I) the value of an 
optimal solution. An optimal fractional solution is a lower bound for the optimal value. We 
denote the optimal fractional solution by LIN {I). 

4.2 Rounding 

We use a rounding technique based on the offline APTAS by Fernandez de La Vega & Lueker 
|13| . As we plan to modify the rounding through the dynamic rounding algorithm we give a 
more abstract approach on how we can round the items to obtain an approximate packing. 
At first we divide the set of items into small ones and large ones. An item i is called small 
if s(i) < e/2, otherwise it is called large. Instance I is partitioned accordingly into the large 
items II and the small items Is- We treat small items and large items differently. Small items 
can be packed using a greedy algorithm and large items need to be rounded using a rounding 
function. We define a rounding function as a function R : II N which maps each large 
item i to a group j. By we denote the set of items being mapped to the same group j, i.e. 
R^ = {i G /l I R{i) = j}- By Aj^ we denote an item i with s{i) = max{s(ifc) | i^ G R^}. Given 
an instance / and a rounding function R, we define the rounded instance I^ by rounding the 
size of every large item i G R^ for j > 1 up to the size s(Aj^) of the largest item in its group. 
Items in R'^ are excluded from instance I^. We write sji{i) for the rounded size of item i in I^. 
Depending on constants c and d, we define the following properties for a rounding function R. 

(A) max{i?(i) I i G II} = c/e"^ for a constant c G 

(B) = \R^\ for all i,j > 2 

(C) \R°\=d\R^\ for a constant d G M+ with d>l 
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(D) s{i) < s{j) ^ R{i) > R{j) 

Any rounding function fulfilling property (A) has at most 0(l/e^) different item sizes and 
hence instance can now be solved approximately using the LP relaxation. The resulting LP 
relaxation has 0(l/e^) rows and can be solved approximately with accuracy (1 + 6) using the 
max-min resource sharing [15 in polynomial time. Based on the fractional solution we obtain 
an integral solution y of the LP with ||y||-^ < (1 + 5)LIN{I^) + C for some additive term C > 0. 
We say a packing B corresponds to a rounding R and solution y if items in R^, . . . , i?™" are 
packed by B according to the integral solution y of the LP. The LP is defined by instance I^. 
Items in R^ are each packed in separate bins. 

Lemma 10. Given instance I with items greater than e/2 and a rounding function R fulfilling 
properties (A) to (D), then OPT{I^) < OPT{I) and < ^eOPT{I). Let y be an integral 
solution of the LP for instance with \\y\\i < (1 + S)LIN{I^) + C for some value C > 0, let 
B be a packing of I which corresponds to R and y and let e' = ^e. Then 

max{Bt{i)} = \\y\\^ + \R°\ < (1 + e' + 5)0PT{I) + C. 

i 

Proof. Let m = max{i?(i) | i G 1^}. Let R^ be the set of items in rounding group i, which 
corresponds to their rounded sizes and let TV' be the set of items in i?*, which corresponds to 
their actual size. Instance contains every item from R^ to i?™, while items from are 
excluded. By property (D) we know, that items in i?* are larger or equal than items in i?*"*"^. 
By property (C) we find for every item in 7i} an unique item in FiP with larger or equal size, 
since the largest item in F{? to which all items are being rounded up is smaller than any item in 
R^. Using property (B) for each item in 7^*+^ we find a unique larger item in R^. Therefore we 
have for every item in the rounded instance L^ an item with larger size in instance / and hence 

OPT{I^) < OPT{I). 

Since the packing B corresponds to a solution y, B gives a solution with maxi {Bt{i)} < 
(1 + 5)LIN{I^) + C + \R^\ bins and since LIN{L^) < OPT{I^) < OPT {I) we obtain that 
maxj{i?((z)} < (1 + 6)0PT(I) + C + l-R*^]. Further, we can bound |-R''|. Since every item in 
/ is of size at least e/2 there is a lower bound for the optimum: OPT{I) > e/2X]o<i<m 1-^*1 — 

e/2Eo<i<m \R^\/d = '^""Vi'^"' ^ ©• Resolving this inequality, we get \R^\ < ^OPT{L) and 
hence \R^\ < e'OPT{I). Since c and d are constant we know \R^\ < e'OPT{I) = 0{eOPT{I)) 
which implies together with the inequality maxj < (1 + 5)0PT{I) + C + l-R^^I < (1 + 

e' + 5) OPT{L) + C. □ 

How can we handle the small items? Actually, small items do not make problems at all. 
We can pack them via FirstFit [5^ on top of the existing large items and still obtain a good 
solution. FirstFit is a greedy algorithm which simply places the current item into the first bin 
having enough space. FirstFit opens a new bin of the item does not fit into any used bin. 

Lemma 11. Let I be an instance with small and large items and given a packing B of the 
large items with maxi {Bt{i)} < K for some K >\. Packing the small items via FirstFit on 
top of packing B gives a new packing of instance I which uses 

max {K,{l + e) OPT {I) + 1} 

bins. 

Given instance I = {ii, . . . ,it}, we define m by m = [1/e^] if [1/e^l is even and otherwise 
m = [1/e^] +1. By definition m is always even. For every instance / we find a rounding function 
R with rounding groups . . . i?™ which fulfills properties (A)-(D) such that \R°\ < 2\R^\ 

and > \R^\. 
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Algorithm 6. 



1. Partition the large items according to the rounding function R in groups ^ • • • 5 R^ 

2. Round up the size of each large item i ^ R^, . . . ,R"^ to s(A^) to obtain instance 

3. Compute a fractional solution x of the LP defined by approximately with ratio (1 + 5) 

4. Round up each component of the fractional solution to obtain an integral solution y for 
the LP for instance 

5. Pack items in R^, . . . , R^ according to the integral solution y 

6. Open a bin for each item i with R{i) = 

7. Pack the small items in Is via FirstFit 

A solution X of with ratio (1 + ^) having m + 1 non-zero components can be computed 



using max-min resource sharing [15 . According to Lemma 10 and 11, the algorithm described 

2d, 
c 



above produces a solution with approximation < {1 + e' + 6) OPT + m + 1 with e' < < 4e. 



4.3 Online Bin Packing 

Let us consider the case where items arrive online. As new items arrive we are allowed to 
repack several items but we intend to keep the migration factor as small as possible. We present 
operations that modify the current rounding Rt and packing Bt to give a solution for the new 
instance. The given operations worsen the approximation but by applying the results from the 
previous section we can maintain an approximation ratio that depends on e. The presented 
rounding technique is similar to the one used in [TD]. In our algorithm we use approximate 
solutions of ILPs in contrast to the APTAS of Epstein & Levin which solve the ILPs optimally. 
Handling with approximate ILPs results in a different analysis of the algorithm because many 
helpful properties of optimal solution are getting lost. 

Note that in an online scenario of bin packing where large and small items arrive online, 
small items do not need to be considered. We use the same techniques as in P^U] to pack small 
items. As a small item arrives we place it via FirstFit [5]. In this case FirstFit increases the 
number of bins being used by at most 1 ([13]) and the migration factor is zero as we repack no 
item. Whenever a new large item arrives several small items might also need to be replaced. 
Every small item in a bin that is repacked by the algorithm, is replaced via FirstFit. Packing 
small items with this strategy does not increase the number of bins that need to be repacked as 
a large item arrives. Later on the migration factor will solely be determined by the number of 
bins that are being repacked. More precisely, we will prove that the number of bins, that need to 
be repacked is bounded by 0{l/e^). Therefore we assume without loss of generality that every 
arriving item is large, i.e. has a size > e/2 (see also [lOj). Our rounding Rt will be constructed 
by three different operations, called the insertion, creation and union operation. The insertion 
operation is performed whenever a large item arrives. This operation is followed by a creation 
or an union operation depending on the phase the algorithm is in. Let I = {ii, . . . ,it} he the 
existing instance as defined above, let R be the corresponding rounding function, let x be a 
fractional solution of the LP generated for the rounded instance I^ and let B be the current 
packing of items in /. We define two subgroups of R^ denoted by R^'^ and R^'^ in the creation 
phase, which are also being modified by the operations. Let /' = IU{ij+i} be the new instance. 
We use the following operations that modify the current rounding R, the packing B and the 
fractional and integral LP solution x and y. We denote with R' , B', x' and y' the new rounding, 
packing and fractional/integral LP solutions for instance I' 
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Ri 



R" 




Insertion Step 

Find the largest j with s(Aj^) > s{it+i). Set R'{it+i) = j and B'{it+i) = B{\f). For every 
= 1, . . . , j we define i?'(Af ) = - 1 and B'{X^) = B{X§_^). Set x' = x and y' = y. 

Modified Insertion Step 

During the creation phase, the algorithm uses the modified insertion operation. Find the largest 
j {j = 1.5 and j = 2.5 included) with s{Xf) > s{it+i). Set R'{it+i) = j and B'{it+i) = B{Xf). 
For every k = l,4,5,...,j we define R'{X^) = k — 1 and B'{X^) = B{X^_^). For every 
= 1.5, 2, 2.5, 3 we define B!{X^) = A; - 0.5 and B'{X^) = B{X^_f^ ^). Set x' = x and y' = y. 

Creation Phase 

The creation phase consists of k creation steps, where k = \R^\. At the end of each creation 
phase we intend to have new rounding groups R^ and R^ created from the subgroups of R^ 
named R^'^ and R^'^. At the beginning of the creation phase we always have \R^\ = 2k and 
R^'^ and R^'^ are empty. In the first step we change the rounding group for all items i with 
-^(^) = J ^ 1 to R'{i) = J ' + 2. Furthermore we say the k largest items of R^ belong to R^'^ 
and the A; smallest items belong to i?^-^. In each of the k creation steps we change the rounding 
function for the largest items A^^g and A^^g. Set R'{X^c^) = 1 and R'{X2,^) = 2. Since items Xf^ 
and A2^5 are moved from R^ to R^ and R^ they have to be covered by the LP. Therefore we 
increase the value of the LP solution by = Xj + 1, x'j = Xj + 1 and y- = + 1, y'j = yj + 1, 
where i,j are defined such that Cj = {1 : ^^/(Afg)} and Cj = {1 : ^/^/(Ag^g)}. For k ^ i,j set 
x'k = Xk and y'^ = yk- 
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R" 



Figure 3: Union operation 

Union Phase 

The union phase consists of k union steps, where k = \B}-\. At the end of each union phase 
we have made out of 4 roundings groups 2 rounding groups with size doubled. For the first 
union step we determine the largest index j with \R^\ < \R^~^^\. If there is no such index then 
set j = m. In each step now set R'{\f ) = j — 1 and R'{\^_2) = j — 3 and for the other 
items i we define R'{i) = R{i). Modify the packing for \^ and Aj^2 by -B'(Aj^) = B{\^_2) and 
place Aj^2 ™to a new bin. Modifying the packing this way implies that we have to change one 
configuration of the fractional and integral LP solution x and y and add one configurations for 
the additional bin. Let Cj be the configuration used by B\Xj). Configuration Cj is replaced 
by a configuration Ci where an item of size SRi{\f_2) is exchanged by an item of size SRi{\f). 
Furthermore we add another configuration Q with an item of size SRi{\f_-^. 

Note that each repacking that we perform in the operations is valid because wc always 
replace items by smaller ones. New packings B' are created in a way that they correspond to 
new integer solution y' . We have to prove that this solutions y' is feasible. Note also, that in a 
creation operation and in a union operation two additional non-zero components of size 1 might 
be created. 

Lemma 12. Applying any operation above on a rounding R and ILP solution y with corre- 
sponding packing B defines a new rounding R' and a new integral solution y' . Solution y' is a 
feasible solution of the LP for instance . 

Proof. We have to analyze how the LP for instance I^' changes in comparison to the LP for 
instance 

Insertion Operation: The right hand side of the LP derived from R' does not change at all 
since the right hand side is determined by the cardinalities \R'^\ = . . . , = \RJ^\. For 

some J > 1 let R^ be the the rounding group where the new item is inserted. By construction 
of the insertion operation for each rounding group R^ with i = 1, . . . there is one item that 
is inserted into group R'^ and one item that is shifted out. Let if be the second largest item of 
rounding group R^. Since the largest item \f in group R^ is shifted to the next group, the size 
sjir{i) of item i in a group R^ is defined by SR'(i) = if. Therefore each item in is rounded 
to the previous smaller value since 3(1^) < s{Xf). Hence configurations of the LP solution for 

can be transformed into feasible configurations for i.e. \\y'\\i = \\y\\i- 
Creation Operation: Note that the rounding groups i?^, for ^ = 1, . . . , m+2 remain identical; 
i.e. R'^ = R^. The groups and get both a new item, but of smaller size. Therefore the 
sizes Sr{i) of all items i ^ are not modified by a creation operation. We have = SRi{i) 

for items in groups R^, . . . , i?™"*'^. Therefore the matrix A = (a((i, j'))) remains the same. Only 
the right hand side h' of the LP from instance is modified (i.e. — = 2). As two new 
configurations arc being added to x and y they cover exactly the enhanced right hand side and 
are therefore a feasible solution of the LP from instance . 

Union Operation: In the union operation we basically change only 4 rounding groups. Sup- 
pose we merge rounding group R^"^ with R>~'^ and rounding group R>-^ with K>. While the 
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size of \R'^~'^\ = \R^~'^\ + 1 and \R'^~^\ = \R^~^\ + 1 is incremented the size of l-R'-'l = \R^\ — 1 
and l-R'-'"^! = — 1 is reduced. Similar to the creation operation, this leads to a change 

in the right hand side of the LP. Two components of the right hand side, which correspond 
to s{X^ ) and s{\^_2) are reduced by 1 and two other components, which correspond to the 
s{\fLi) and s{X^L^) are increased by 1. Furthermore the sizes of items in R'^ and R'^^"^ are 
equal or smaller than the sizes of items in R^ and R^~'^ since s_R/(i) = if for all items i G R'^ 
and sji'{i) = if_2 for all items i S R!^~'^. Xf and Aj^2 shifted to the next rounding groups. 
Consider a feasible configuration C of the LP for instance I^. Then the modified configuration 
C (with replaced item sizes) is also feasible in the LP for instance . The new solutions x' 
and y' use the modified configurations and cover the right hand side of the LP. □ 

The operations are used as described in Algorithm [7] below. We apply the algorithm on 
a rounding function Rq and instance Iq = {ii, . . At}- We suppose that |i?Q| = \Rq\ = . . . = 
\Rq^\ = K for some K > Q and hence T = K{m + 1). An improve(a, x, y, 5) statement stands 
for a call of Algorithm [5] with improvement a = a, fractional solution x and integral solution 

y. The variable part x'"'^'^ is defined by x'""''^ = 2 ^^V^,^'^^ x. After an improve call the packing is 

ll^lli 

changed according to the new integral solution. Since during a creation operation and a union 
operation two additional non-zero components of size 1 might appear, we change the parameter 
i of Algorithm [5] slightly to I'. Parameter ^' is defined maximally such that the sum of the 
smallest components y'l, . . . , are X]i<i<£ v'ai ^ (m + 2)(l/5 + 2) + 2. The two additional non- 
zero components belong to components yi , . . . yi+2 and are therefore reduced in step 6 along 
with the others. 

Algorithm 7. 



for i := 1 to K do 

get new item; 
improve (1 ,x,y,6); insert; 
for i := 1 to m/2 do 

/* Creation Phase */ 
for j := 1 to K do 

get new item; 

improve(l, X, y, 5 ); 

modified insert; 

create; 

/* Union Phase */ 
for j := 1 to K do 

get new item; 
improve (2, x,y,5); 
insert; 
union; 



In the following we present how the algorithm changes the rounding groups for m = 6. The 
table presents the state of each rounding groups after each phase. One can see that after the 
execution of Algorithm [7] each rounding group has exactly 2K items. We prove the general case 
for arbitrary m: Every rounding has exactly 2K items after the execution of Algorithm [7| 

Lemma 13. Let Rq be the rounding function at the beginning of the algorithm. Suppose that 
every rounding group i?Q,...,i?g^ has exactly K items. Then after the execution of the algo- 
rithm above the computed rounding function R^ after T insertions has m + 1 rounding groups 
R%,...,Rf with \R^j,\ =2Kfore = 0,...,m. 
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Proof. The algorithm starts with a rounding function that contains exactly T items. After the 
first K insertion steps rounding function Rk is of the form: \R^j^\ = 2K, \ R\^\ = K, . . . , |-R^| = 
K since K items are shifted to R^ while the cardinalities of the other rounding groups remain 
the same. During the next K arrivals, the algorithm is in the creation phase. We perform a 
creation operation after each insertion. For each item shifted to R^, two items are shifted to the 
new created groups R^ and R^. At the end of the first creation phase, the rounding function 
R2K satisfies = K,\RIj^\ = K,...,\R^+^\ = K. In the following union phase, the 

rounding groups i?™^'^, ij^+i and i?"*, R"^~^ are merged together. For each union operation, 
one item is shifted from i?™"*"^ to i?™"^^ and another from i?™" to R"^~^. Since there are K 
insert operations in the union phase, rounding group \R^k\ ~ ^-^J-^sicI = K, . . . ,\R^~'^\ = 
K, \R]^~^\ = 2K, = 2K. After the next creation and union phase, the number of rounding 
groups is also m + 1. On the other hand we have two additional groups l-Rs^^^l = l-^^^l = 2-?^. 
After j < m/2 creation and union phases the rounding function R2jK+K is by induction of the 

frirm \ ff^ \ — 9 K \ \ — K \ r)m-2j i _ ^ |„m-2j+l| _ n jy- \ nm \ — n ly- 

lOim \-n2jK+K\ — \-"-2jK+K\ —-"■)■■■) \-'^2jK+K\ ~ \^2jK+K I — ^-tv, • • ■ , \n2jK+K\ — 

This can be proved by induction on j. For j — m/2 — 1 we get |-R2m_ft'— E"! — ^-ftT, |-R2mi<'— K I — 
K,\R2mK-K\ ~ ^-•\^2mK-K\ ~ "^^^ ' ' ■ \^2mK-K\ ~ ■ Therefore, after one additional 
creation and union phase, we obtain m + 1 groups \R2mK+K\ ~ for £ = 0, . . . , m + 1. □ 

Using Algorithm [7] with a starting rounding function Rt that has m rounding groups and the 
property that \R^\ = \Rt\ = . . . = |-Rj^| produces according to the lemma a rounding function 
i?2r that has also m rounding groups of equal size, but with cardinality doubled. Therefore 
we can use Algorithm [7] repetitively to always get suitable rounding functions. The following 
algorithm is our final online AFPTAS for the classical bin packing problem. Let St be the sum 
of all item sizes of items ii, . . .it. 

Algorithm 8. 

• While St < (m + 2) (1/5 + 4) and {m + 1) does not divide t get the new item i^+i and use 
the offline AFPTASlm with an LP of approximation ratio (1 + 5). 

• Afterwards use Algorithm^ repetitively to obtain a packing for each instance 

By using the offline AFPTAS for small instances we can make sure that Algorithm [7] is 
started with a suitable rounding function. Since Algorithm [7] always produces a rounding 
function fulfilling properties (A) to (D) and m + 1 divides the current number of items t, every 
rounding group R^,...,R^ has the same number of item sizes as the algorithm leaves the 
while-loop in the first step. 

In the following we give a bound for the rounding functions Rt that we produce in every 
step of the algorithm. It remains to prove that the approximation during the execution of 
Algorithm [7] can be bounded. Therefore we define a relation between rounding functions. Let 
R and R be two rounding functions, with R having rh rounding groups for some fh G ©(l/e^). 
We can embed R into R in symbols R < R, if \R^\ < \R^\ and for every item i G I\R^ we have 
< Sjiii). A relation R<R always implies that + OPT{I^) < + OPT{I^). 
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Lemma 14. For each t G N^, we can embed Rt into a function Rt, which fulfills properties (A ) 
to (D). Rounding function Rt has parameter c > 1/4 for property (A) and d < 2 for property 
(C). 

Proof. Since we basically shift largest items to the following rounding group we designed op- 
erations insertion, creation and union in a way that property (D) is never being violated by 



any Rj for j < t. As shown in the proof of Lemma 12 the number of rounding groups remains 
constant between m + 1 at the end of a union phase and m + 3 during the creation and union 
phase. Suppose the algorithm above has started with a number of items T < t < 2T and 
rounding groups Rj^, . . . i?™, which are being modified by the algorithm. We define the rounding 
function R = Rt in which Rt can be embeded in the following way: Function R has rounding 
groups . .,fiL2i?J-i with = . .. = |^L27fJ-i| = 2K and |^°| = 2K + {t mod 2K). Since 
every rounding group of R except R^ has the same number of items, the rounding function R 
fulfills property (B). Rounding function R fulfills property (C) because 2K < \R^\ < AK and 
l^^l = 2K. This implies constant d <2. We prove property (A) by giving an upper and a lower 
bound for maxj that are both in Q{^). Recall that T / K = m + l. On the one hand we 

get maxi [R{i)] = [^\ - I < ^ - I < ^ - I < {m + I) - I = m < \^'\ + I < ^ + 2. On 

the other hand max^ = [^\ - 1 > [^J - 1 > L*^] -1>t-1 = 2??-1- ^ince 

e < 1/2 weg get c > 1/4. It remains to prove that we can embed Rt in R i.e. Rt < R- Since R^^ 
never exceeds 2K items and 2K < AK we get < \R^\. According to the proof of Lemma 12 



and the construction of the creation operation, Rt is during the creation phase of the following 
form: = 2K - a, \Rl\ = a, \Rl\ = a, = . . . = \r{\ = K, = . . . = 2K for 

some a < K and j < m + 2. Rounding function R has in every rounding group R^ for j > 1 
exactly 2K items. Since property (D) holds for both rounding function Rt and i?, the rounding 
groups . . . , contain the same items as the last m + 2— j rounding groups of R. Items 

in these groups are therefore rounded identically. For some fh let i?™ be the rounding group 
which contains the same items as R{^^ ■ Since rounding groups i?f , . . . , i?j each contain exactly 
K items, the rounding groups i?^, . . . , contain the items of exactly two rounding groups. 

Therefore the items in . . . , R^~^ are rounded to a smaller size compared to using R. Items 
that belong to R^ and -R^ are contained in R^ and by definition do not need to be considered. 
Hence, any rounding function Rt which is in an creation phase can be embedded into an R. By 



construction of the union operation and the proof of Lemma 12 -Rt is during the union phase 



of the form \R^^ \ = K + a, \Rl\ = ... = \Ri~^\ = K, \r{~^\ = K + a, \r{~'^\ = K - a, \Ri~^\ = 
K + a, \R{\ = K-a, \Ri^^\ = ... \RT~^'^\ = for some a< K and j < m + 2. As shown in the 
union phase, items in Rt~^^ , • • . , R^'^'^ are rounded equally in R. As the sum of i?r and Rl IS 
2K and the sum of R{~^ and R^""^ is 2K the items of R{~^ and R^ and the items in R{~^ and 
belong in R to the same rounding group and are hence rounded equally or to a smaller 
size compared to using R. Items in , . . . , R?t~^ are each of size K and are rounded equally or 
to a smaller size than using R since the same argument as in the creation phase holds. □ 



Define e by e = ^e. As Rt fulfills property (A) to (D), we obtain by Lemma 10 and Rt < R 



the following two equations for every t: 

1. OPT{I^') < OPT {It) 

2. \R^^\ < \R^^ \ < ^eOPT{It) < eOPT{It) 

Note that since e < 1/2 we have e < ^. Recall that Rt < Rt implies that \R^\ + OPT{I^^) < 
\R^\ + OPT{I^^) and that LIN{I^^)+m > OPT{I^') (rounding up a basic feasible solution). Let 
us discuss how the methods from the previous section apply to the presented online algorithm. 
The procedure improve is implemented by using Algorithm [5] in order to get an improved 
solution for instance Algorithm [5] is applied using 6 as the approximation parameter. In 
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the following lemma we prove that applying Algorithm [5] to improve a solution for impacts 
the overall approximation A = e + 5 + e6 in the same way. We define C = AOPT(It) + m. 

Theorem 15. Given a rounding function Rt and an LP defined for Let x he a fractional 
solution of the LP with \\x\\^ + \R^^ \ < {l + A)OPT{It) and \\x\\-^ > 2a(l/l+l) and \\x\\^ = (1 + 
5')LLN{I^^) for some 5' > 0. Let y be an integral solution of the LP with \\y\\i > (m+ 2) (1/5+2) 
and corresponding packing B such that maxji?j(i) = ||y||^ + \R^\ < (1 + 2A)0PT{Lt) + m. 
Suppose x and y have the same number < C of non-zero components and for all components i 
we have yi > Xi. Then using Algorithm^on x andy returns new solutions x' with + < 
(1 + A) OPT(Lt) — a and integral solution y' with corresponding packing B[ such that 

max5^(i) < {l + 2A)0PT{It) +m- a. 

i 

Further, both solutions x' and y' have the same number < C of non-zero components and for 
each component we have x[ < y[. 

Proof. As shown in the following, Algorithm [5] maintains the property that x and y have the 
same number of non-zero components and that Xi < yi since we can use Theorem[8]and Corollary 
|9| By condition we have maxj Bt{i) = \\y\\^ + \R^\ < (1 + 2A)0PT{It) + m. Since OPT{Lt) < 
OPT{I^) + \R^^ \ we obtain for the integral solution y that ||y||^ < 2A0PT{It)+m+ OPT{I^) < 
2A0PT{It) +m + LIN{I^) + m. Hence by definition of C we get ||y||^ < LIN{L^) + 2C . This 
is one requirement to use Theorem [8] or Corollary [9| We look at the cases separately where on 
the one hand 5' < 6 and on the other hand 5' > 5. 

Case 1, (5' < (5: At first we give an upper bound for LIN{I^^): We get LIN{I^-^) < 
OPT{I^') < OPT{I^') < OPT{L^^) + \Rt^\ < {l + e)OPT{It) using that Rt < Rt- This 

implies that 5LIN{I^*) < 50PT{L^^) < {5+5e)OPT{It) < C. Algorithm|5]returns by Theorem 
[sja solution with ||a;'||^ < (1+(5)L/A'^(/^')— a and an integral solution y' with \\y'\\i < ||x'||-^+C 
or lly'lli < For the term + we get + |/??| < {1 + 6)0PtIi^') - a + \R^\. 

Using that Rt can be embedded in Rt we get + OPT{L^') < |i?t°| + OPT{I^') < 
OPTpt) + eOPT{It). Therefore \\x'\\^ + < 60PT{I^') -a + OPT{It) + eOPT{Lt) < 
{5+6e)OPT{It)-a + {l + e)OPT{It) < {l+A)OPT{Lt)-a. In the case where \\y'\\^ < \\x'\\-^ + C 
we can bound the number of bins of the new packing B' by maxji?^(i) = + |^t| < 

ll^^'lli + \Rt \ + C < (1 + A)OPT{It) - a + C = (1 + 2A)0PT{Lt) + m - a. In the case that 
lly'lli < ||y||i — a we obtain maxji?j(i) = \\y'\\i + \R^\ < \\y\\i — a + \Rt \ = maxji?j(z) — a < 
{l + 2A)0PT{Lt)_+m-a. 

Case 2, 6' > 5: By condition we have \\x\\-^ + \R^\ < (1 + A)OPT{Lt). Since OPT{Lt) < 
OPT{I^') + \R^\ we obtain for the solution X that \\x\\-^ < AOPT{Lt) + OPT{I^') < AOPT{It) + 
LIN{L^^) + m. Hence by definition of C this implies ||x||i < LIN{I^^) + C and therefore 
6' LIN{I^*) < C, which fulfills the requirements of Corollaryjoj Using Algorithm [H] on solutions 
X with ||x||i = (1 + 5')LIN{L^^) and y with ||y||^ < LIN{I^) + 2C we obtain by Corollary [o] 
a fractional solution x' with \\x'\\-^ < \\x\\-^ — a and an integral solution y' with either \\y'\\i < 
llyll;^ — a or < ||x||]^ + C — a. So for the new packing 5' we can guarantee, that maxj = 
+ < ||y||i-a + |i?j°| = maxi Bt{i) - a < {l + 2A)0PT{Lt) + m- a if \\y'\\^ < ||y||i-a. 
If lly'lli < + C — Q, we can guarantee that maxj B^{i) = \\y'\\i + \Rt \ < II^^Hi + | +C — a < 
(1 + A)OPT{It) + C - Q < (1 + 2A)0PT{It) + m-a. Furthermore we know by Corollary |9] 
that x' and y' have at most C non-zero components. □ 

Set (5 = e. Then A = 2e + = ©(e). We get the central theorem: 

Theorem 16. Algorithm^is a fully robust AFPTAS for the bin packing problem. 

Proof. While instances are small Algorithm |8] uses the offline AFPTAS (see Algorithm [g]). 
Using Algorithm [g] we get a packing Bt for instance It that uses at most maxj Bt{i) < (1 + e' + 
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5) OPT{It) + ^ + 1 bins, where e' < 4e < e. Since the instance is small the migration factor 
is bounded although we might repack every single item. Let r be the first index where the 
algorithm leaves the while-loop. By condition we are in the while loop while St < (m-|-2) (1/5+4) 
and t does not divide m + 1. Hence Sr < {rn + 2) (1/^ + 4) + m = 0{l/e^). The migration 
factor for instances It with t < t \s therefore bounded by "^St = 0{l/e^) since every arriving 
item has size at least e/2. The approximation guarantee for small instances is bounded by 
maxj Bt{i) < {1 + 5 + e) OPT{It) + m + 1. In the following we consider large instances It with 
t > T. 

Full robustness: The migration factor for some consecutive packings Bt and Bt+i is 
bounded by the migration of the improve-call plus the migration of an insertion and an union 
operation. The operations create requires no shifting of items at all. As proven in the previous 
section, an improve-call changes at most 0{m/5) components of a solution y. Since the arriving 
item is large with size > e/2, changing a complete configuration requires migration of at most 
0{l/e). Combined this results in a migration factor for the improve-call 0(m/A^) = 0{l/e^) if 
we use Algorithm |5j By construction of the insertion operation it shifts in worst case one item 
per rounding group. Having ©(l/e^) rounding groups this gives a migration factor of at most 
0(1/6^). Therefore the complete migration is bounded by 0{l/e^). 

Running time: The running time is dominated by the max-min resource sharing (see 
Algorithm [5]) and the number of non-zero components. The number of non-zero components 
is bounded by AOPT{It) + m < At + ^ + 1 and is therefore polynomial in ^ and t. As the 
running time for the max-min resource sharing is also polynomial in ^ (see [15])) the running 
time is clearly polynomial in t and j. 

Approximation: We prove by induction that four properties hold for any packing Bt and 
corresponding LP solutions. Given fractional solutions x and integral solution y of the LP 



defined by instance /^*. Properties (2)-(4) are necessary to apply Theorem 16 and property (1) 
provides the wished approximation ratio for the bin packing problem. 

1. packing Bt uses at most (1 + 2A) OPT(It) + m bins 

2. \\x\\^ + \R^\ < {l + A)OPT{It) 

3. for every configuration i we have Xi < yi 

4. X and y have the same number of non-zero components and that number is bounded by 
AOPT{It)+m 



To apply Theorem 15 we furthermore need a guaranteed minimal size for ||x||^ and \\y\ 



1- 

According to Theorem 15 integral solution y needs ||y||^ > (m + 2) (1/6 + 2) and > 

4(1/(5 + 1) as we set at most a = 2. By condition of the while-loop we know that any in- 
stance St > {m + 2){l/6 + 6)^ Since OPT{It) < \\y\\-^ + < |[y||i + eOPT{It) we get 
||y|li > (1 -e)OPT{It) = (1 - 5)0PT{It).By OPT{It) > (m + 2)(l/(5 + 4) we finally get that 
||y|li > (l-5)(m+2)(l/5+6) > (m+2)(l/5+6)-(m+2)(l + 65) > (m+2)(l/5+6)-4(m+2) = 
(m + 2){l/6 + 2). Since OPT{It) < \\x\\^ + m + \R^\ we obtain by the same argument that 

Iklli > (m + 2){l/6 + 2) - m > (m + 2){l/5 + 1) and since m = l/e > 1/5 > 2 we get that 
Iklli > 4(1/5 + 1). 

In the case that t = t we have by the offline algorithm that the number of non-zero com- 
ponents = m + 1 < AOPT(It) + m since OPT {It) > St > 1/A. The number of used bins is 
bounded by max^ Bt{i) < (1 + 5 + e) OPT{It) + m + 1 < (1 + 2A) OPT{It) + m (note e' < e) 
and property (2) is fulfilled for the same reason. Furthermore in the offline algorithm every 
component Xi is rounded up to obtain the integral component y^. Therefore all properties (1)- 
(4) are fulfilled for t < t and the induction basis holds. Now let Bt be a packing for t > t for 
instance It with solutions x and y of the LP defined by I^* . Suppose by induction that property 
(l)-(4) hold. We have to prove that these properties also hold for Bt+i and the corresponding 
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solutions of the LP defined by 1^^+^. Packing -Bt+i is created by using an improve call for x 
and y followed by an insertion operation and optional, an union or a creation operation, 
improve: Let x' be the resulting fractional solution of Algorithm |5] let y' be the resulting 
integral solution of Algorithm |5] and let B'^ be the corresponding packing. Properties (l)-(4) 
are fulfilled for x, y and Bt by induction hypothesis. Hence we can use Theorem [T5| By 



Theorem 15 properties (l)-(4) are then still fulfilled for x', y' and B[ and moreover we get 
Ik'lli + |^< (1 + A)OPT{It) - a and maxiS^(i) < (1 + 2A)0PT{It) + m - a for a = 2 or 
a = 2. 

operations: First we take a look at how the operations modify \\y'\\i and \Rf\. By 

construction of the insertion operation, the LP solutions x' and y' are not modified while 
increases by 1. By construction of the creation operation \\x'\\i and are increased by 2 and 
l-R^I decreases by 2. By construction of the union operation, \\x'\\^ and ||y'||^ are increased by 1 
and I remains constant. Property (1): Let x" be the fractional solution and y" be the integral 
solution after using operations on x' and y' . Packing Bt+i equals maxj B^+i = \\y"\\i + l-R^+il- 
According to the operations an insertion operation yields maxj = + + 1 = 

maxj i?^ + 1. An insertion operation followed by an union operation yields maxj B^+i = 
||y'||-i^ + 1 + l-R^I + 1 = maxj i?j + 2 and an insertion operation followed by a creation oper- 
ation yields maxj Bt+i = ||y'||]^-|-2-|-|ii^| — 1 = maxj B[-\-l. Algorithm [t] is designed that in the 
union phase maxj B'^ < {1 + 2A) OPT(It) +m — 2 since there is an improve call with a = 2 and 
otherwise maxj B[ < (1 + 2A) OPT{It) + m — 1 since there is an improve call with a = 1. There- 
fore we have in any case that .Bt+i uses at most {\ + 2/\)0PT{It)+m < {l + 2A)0PT{It+i) + m 
bins. The proof that property (2) holds is symmetric since \\x'\\i increases in the same way as 
lly'lli and + \R^\ < {I + A)OPT{It) - a for a = 1 or a = 2. For property (3) note that 

in the operations a configuration Xi of the fractional solution is increased by 1 if and only if a 
configuration yi is increased by 1. Therefore the property that for all configurations < y'.' 
retains from x' and y' . By Theorem 15 the number of non-zero components of x' and y' is 
bounded by AOPT(It) +m < AOPT(It+i) +m. By construction of the creation operation and 
union operation x" and y" might have two additional non-zero components. But since these 
are being reduced by Algorithm [s] (note that we increased the number of components being 
reduced in step 6 by 2), the LP solutions x" and y" have at most AOPT{It-^-i) + m non-zero 
components which proves property (4). □ 



4.4 Running Time 

Storing items that are in the same rounding group in a heap structure, we can perform each 
operation (insertion, creation and union) in time 0{-^ log(e^t)). Furthermore Algorithmjs needs 
to look through all non-zero components. The number of non-zero components is bounded by 
O{e0PT) = 0{et). Main part of the complexity lies in finding an approximate LP solution. 
Let M{n) be the time to solve a system of n linear equations. The running time of max-min 
resource sharing is then in our case 0{M{^)^ + 4) (see ^). Therefore the running time of 
the Algorithm is 0{M{^)^ + et + ^ \og{eH)). 



5 Conclusion 

Based on approximate solutions, we developed an analogon to a theorem of Cook et al. [7^. Our 
improvement helps to develop online algorithms with a migration factor that is bounded by a 
polynomial in 1/e, while algorithms based on Cook's theorem usually have exponential migration 
factors. We therefore applied our techniques to the famous online bin packing problem. This led 
to the creation of the first fully robust AFPTAS for an NP-hard online optimization problem. 
The migration factor of our algorithm is of size O(^), which is a notable reduction compared 
to previous robust algorithms. When a new item arrives at time t the algorithm needs running 
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time of 0{M{-^)-^ + et + ^log(e^t)), where M(n) is the time to solve a system of n Hnear 
equations. Any improvement to the max-min resource sharing algorithm based on the special 
structure of bin packing would immediately speed up our online algorithm. We believe that 
there is room to reduce the running time and the migration factor. Note for example that we 
give only a very rough bound for the migration factor as the algorithm repacks 0{^) bins. 
Repacking these bins in a more carefully way might lead to a smaller migration factor. An open 
question is the existence of an AFPTAS with a constant migration factor that is independent of 
e. We mention in closing that the LP /ILP-techniques presented are very general and hence can 
possibly be used to obtain fully robust algorithms for several other online optimization problems 
as well (i.e. multi-commodity flow, strip packing, scheduling with malleable/moldable task or 
scheduling with resource constraints) . 
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